Draft from February 2, 2008, to appear in ApJ 



High-spatial-resolution CN and CS observation of two regions of 

massive star formation 

H. Beuther 

^ ! Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 

(N 

hbeutherOcf a . harvard . edu 
P. Schilke, F. Wyrowski 
Max-Planck-Institute for Radioastronomy, Auf dem Hiigel 69, 53121 Bonn, Germany 



0^ 



O 



> 

^ : ABSTRACT 

m 

O ! Molecular line CN, CS and mm continuum observations of two intermediate- 

to high-mass star-forming regions - IRAS 20293+3952 and IRAS 19410+2336 - 
obtained with the Plateau de Bure Interferometer at high spatial resolution reveal 
Qh| interesting characteristics of the gas and dust emission. In spite of the expectation 

O \ that the CN and CS morphology might closely follow the dense gas traced by 

^ \ the dust continuum, both molecules avoid the most central cores. Comparing 

the relative line strengths of various CN hyperfine components, this appears not 
to be an opacity effect but to be due to chemical and physical effects. The CN 
^ ' data also indicate enhanced emission toward the different molecular outflows in 

d \ the region. Regarding CS, avoiding the central cores can be due to high optical 

depth, but the data also show that the CS emission is nearly always associated 
with the outflows of the region. Therefore, neither CS nor CN appear well suited 
for dense gas and disk studies in these two sources, and we recommend the use of 
different molecules for future massive disk studies. An analysis of the 1 and 3 mm 
continuum fluxes toward IRAS 20293+3952 reveals that the dust opacity index 
j3 is lower than the canonical value of 2. Tentatively, we identify a decreasing 
gradient of [3 from the edge of the core to the core center. This could be due to 
increasing optical depth toward the core center and/or grain growth within the 
densest cores and potential central disks. We detect 3 mm continuum emission 
toward the collimated outflow emanating from IRAS 20293+3952. The spectral 
index of a ~ 0.8 in this region is consistent with standard models for collimated 
ionized winds. 
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Subject headings: star: formation - techniques: interferometric - astrochemistry 
- ISM: individual (IRAS 20293+395; IRAS 19217+2336) 

1. Introduction 

Massive star-forming regions consist of hundreds to thousands of Mq of dense gas from 
which future clusters are forming. The local densities and chemistries vary significantly 
throughout the regions, and if one wants to observe the densest regions, and in particular 
potential disks at the core center, the right molecular transitions need to be chosen to study 
the gas in detail. Especially the question of disks is still open in spite of large efforts 
over recent years. Pure continuum observations lack the necessary velocity information 
to prove real disk emission, maser emission often only allows ambiguous interpretations, 
and thermal line emission of dense gas tracers is often contaminated by emission from the 
ambient gas (e.g., Cesaroni et al. 1997; Norris et al. 1998; Walsh et al. 1998; Zhang et al. 
2002; Sandell et al. 2003). Nevertheless, we think that observing thermal emission from the 
right molecules - which trace disk-like structures and are weak in ambient gas emission - 
is the most promising path to observe and understand potential genuine disk emission in 
high-mass star formation (e.g., Cesaroni et al. 1997; Zhang et al. 2002; Beltran et al. 2004). 
The question remains: which molecular transitions are best suited to pursue this case? 

Many single-dish molecular-line studies of massive star-forming regions have been un- 
dertaken (e.g., Schilke et al. 1997; Hatchell et al. 1998; Thompson & Macdonald 2003), but 
only very few regions have been observed at the high spatial resolution necessary to disen- 
tangle the densest molecular clumps from the potential disks (Cesaroni et al. 1999; Sandell 
et al. 2003). Here, we present high-spatial-resolution Plateau de Bure Interferometer (PdBP) 
observations of two intermediate to high-mass star-forming regions in two different potential 
dense-gas and disk-tracing molecules, CN and CS. 

Observations toward low-mass T Tauri stars have shown that CN is enhanced in their 
disks (Dutrey et al. 1997), and theoretical work by Aikawa & Herbst (2001) suggested that 
CN could be a good tracer of gaseous disks in deeply embedded sources as well. Therefore, we 
observed the intermediate-mass region IRAS 20293+3952 with the PdBI in several different 
transitions of CN in the 1 and the 3 mm band to test the suitabihty of CN as a potential 
massive disk tracer. The mm sources we are interested in are part of an intermediate- 
mass sub-cluster in the near vicinity of an ultracompact Hll region which provides most of 
the luminosity within the region. For more details on the source and its multiple outflow 
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properties we refer to Beuther et al. (2004b) and Table 1. 

The other PdBI observations aimed at the high-mass star-forming region IRAS 19410+2336 
in two CS transitions and the continuum at 1 and 3 mm. CS has high critical densities 
(> lO^cm"^) and is expected to follow the dust emission well (Shirley et al. 2003). These 
observations were meant to study the large scale CS gas and velocity distribution as well as 
tracing the innermost cores and possible disks. IRAS 19410+2336 has been studied inten- 
sively over recent years, and the complex outflow and mm continuum properties are outlined 
in Beuther et al. (2004b) and Beuther & Schilke (2004). Further source parameters are listed 
in Table 1. 

Section 2 describes the technical details of the observations, and in section 3 we present 
the results obtained for both sources. Finally, section 4 discusses these results in the frame- 
work of future dense gas and disk studies in high-mass star formation. 

2. Observations 

We observed IRAS 20293+3952 and IRAS 19410+2336 in several runs over the two win- 
ters 2001/2002 and 2002/2003 with the Plateau de Bure Interferometer (PdBI) at 1 and 
3 mm. The observing parameters for the various observations are listed in Table 1. 

The typical system temperature at 1.3 mm was about 300 K and at 3 mm about 135 K. 
The phase noise at both wavelengths was mostly below 20° and always below 30°. At- 
mospheric phase correction based on the 1.3 mm total power was applied. For continuum 
measurements, we placed two 320 MHz correlator units in each band to cover the largest pos- 
sible bandwidths. The final continuum images were produced using only the line-free parts 
of the spectrum. We estimate the final fiux density accuracy to be ~ 15%. Synthesized 
beams at the different wavelengths are hsted in Table 1. 

3. Results 

One of the goals of both observations - CN in IRAS 20293+3952 and CS in IRAS 19410+2336 - 
was to trace the densest gas and in particular the massive disks at the very centers of the 
evolving protoclusters. Before going into details, we note that at the highest spatial res- 
olution, in both cases CN and CS do not trace the central condensations, but mainly gas 
in their near vicinity. Thus, the data show that CN and CS are not disk tracers in these 
two sources. This is distinctively different from the observations of low-mass T Tauri disks 
(Dutrey et al. 1997). 
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3.1. Millimeter continuum emission in IRAS 20293+3952 

The three mm peaks (mml to mm3) observed with the PdBI at 2.6 mm (Figs. 1, 2 ) 
correspond well with the three peaks previously observed by Beuther et al. (2004b) at 1.3 mm. 
The only differences are that the previous 1.3 mm data showed a small sub-component north- 
east of the main peak mml, whereas this new data show an elongation to the south-west of it 
(mmla, Fig. 2). The north-eastern 1.3 mm sub-component of the previous dataset could be 
due to dust emission which might be filtered out by the larger configurations used to observe 
the source this time (AB of PdBI compared to CD in the previous observation). Furthermore, 
at 3 mm the sensitivity for dust emission decreases substantially {S oc z/^+^ with /? being the 
dust opacity index), and therefore we simply might not detect this sub-component at the 
lower frequency. 

Contrary to this, the sub-component mmla and its elongation in this new dataset is 
spatially correlated with the collimated outflow in this direction (arrows in Fig. 1, Beuther 
et al. 2004b). This is similar to the recent observations of LI 157 where Gueth et al. (2003) 
distinguish dust emission due to the outflow from the protostellar dust emission. It is likely 
that the sub-component mmla is caused by the collimated outflow. We will discuss the 
scenario of a collimated ionized wind with low spectral index in §4.2. 

Assuming optically thin dust emission for the mm continuum emission, we can estimate 
the dust masses of mml to mm3. Setting the dust opacity index /3 to 1 (see the following 
paragraph) and using a dust temperature of 56 K (Sridharan et al. 2002), the derived core 
masses from the 3 mm data are between 2.8 and 0.4 Mq (Table 2), close to the values derived 
previously by Beuther et al. (2004b). For mml we can also derive a mass from the 1.3 mm 
data which is about a factor 2 lower. This difference can be attributed to the stronger spatial 
filtering effects of the interferometer at shorter wavelengths. 

Spectral indices: The high spatial resolution and good signal-to-noise ratio of the 
continuum data allow us to analyze the spectral index in the mm regime. However, due 
to the different wavelengths, the uv-coverages and thus spatial resolution of the datasets 
are different (Fig. 2). To accurately compare the continuum data it is crucial to employ 
the same similarly sampled uv-range (in kA) in both bands. Therefore, we flagged all data 
outside the uv-range 13 — 150 kA (covered by both frequency bands) and re-imaged both 
datasets with the same beam of 1.5". These new datasets are now well suited to determine 
spectral indices within the central mm core. 

The right panel in Figure 2 presents the spectral index map of the central core mml 
assuming a power law S oc z/° in the Rayleigh- Jeans regime with a = 2 + f3 where f3 is 
the dust opacity index. Since the calibration uncertainty is 15%, we estimate the error in 
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a and /3 to be ~ 0.3. We find that the spectral index a decreases from ~ 3.4 at the outer 
edge of the core to values as low as ~ 2.1 at the very center of mml. Assuming a constant 
temperature throughout the core, this corresponds to dust opacity indices (3 decreasing from 
1.4 at the core edge to 0.1 at the center. The variation of a is more pronounced in north-west 
south-eastern direction perpendicular to the outflow than in the direction of the outflow. It 
is likely that the outflow affects the mm continuum emission of the core in such a way that 
the observed a in the outflow direction is lower than that perpendicular to it (see also the 
end of this section and 4.2). 

Figure 3 outlines a different approach analyzing the data directly in the uv-plane. The 
1 and 3 mm data can both be fitted with a Gaussian plus a point source (excluding the 
data-points < 20 kA) . This way, we derive spectral indices of the Gaussian and point source 
components of 2.7 and 1.6, respectively. Assuming that the Gaussian component is from a 
dust core with constant temperature this results in a dust opacity index (5 of 0.7. The central 
point source could be either a small optically thick disk or an unresolved hypercompact Hll 
region. 

In addition to the main core mml, we can estimate the spectral index a for the outflow 
associated emission feature (Table 2). We flnd a lower spectral index toward mmla with 
Q;(mmla) ~ 0.8. Thus, a is lower in the outflow region than within the central core. 

3.2. CN in IRAS 20293+3952 

Figure 1 presents the CN(l-O) and 2.6 mm continuum emission as contours overlaid on 
the single-dish large-scale 1.2 mm image obtained with the MAMBO-array (Beuther et al. 
2002b). The small offset between the single-dish and the interferometric mm peaks is likely 
due to pointing inaccuracy of the single-dish data estimated to be 5" (Beuther et al. 2002b). 

On large scales, the CN(l-O) emission follows the mm dust emission outlined by the 
MAMBO 1.2 mm data. But at smaller spatial scales we flnd that the CN peaks are always 
offset from the mm peaks. Because the interferometric CN line and mm continuum data are 
taken simultaneously, the relative positional accuracy between the various maps is extremely 
high and the offsets are real. We observe similar morphologies in the 1.3 mm band (Fig. 4). It 
appears that the CN(2-1) emission forms nearly a ring-like structure around the continuum 
source. While in the 1.3 mm band only one line is strong enough to be imaged in these 
observations, in the 2.6 mm band we have the main line plus 6 satellite lines due to the 
hyperflne splitting of CN (e.g., Simon 1997). Figure 5 shows a spectrum of the main line 
and one hyperflne component (integrated over a region of AR.A.[— 5, 1]" and ADec.[— 5, 1]"). 
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The ratio of the peak intensities of the two hnes is of the order of the theoretical value of 
~ 2.7 (Simon 1997). To first order, this indicates that all lines are optically thin and trace 
the whole CN column densities. Therefore, the weak CN emission toward the mm peaks is 
a real physical effect and not due to the optical depth of the observed lines. 

One would also like to derive excitation temperature and column density estimates 
from the CN(2-1) and (1-0) transitions. However, missing flux strongly deteriorating the 
temperature and column density estimates is a problem. Therefore, we refrain from further 
analysis of physical parameters analysis based on this data. Short spacing information is 
crucial to proceed on this path. 

3.3. CS in IRAS 19410+2336 

The mm continuum data and their interpretation are presented in a separate paper 
(Beuther & Schilke 2004). Although we filter out about 80-90% percent of the continuum 
emission, the spatial filtering is even more severe for the line emission. Comparing the 
interferometric measured fiuxes (see example spectra in Figure 6) with the values derived 
from single-dish observations (Beuther et al. 2002b), we find that of the order of only a few 
percent of the total CS fine emission is recovered in the interferometer data. Therefore, the 
dust continuum emission appears to be more strongly associated with the dense cores than 
CS. 

Figure 7 (left panel) presents the integrated PdBI CS(2-1) observations {Av [17,27] 
km/s) as a contours overlaid on the single-dish 1.2 mm dust emission (Beuther et al. 2002b). 
To associate the CS emission with the many outflows in the regions, we show with arrows, 
ellipses and markers the 7 (or even 9?) outflows in the region, as presented in Beuther 
et al. (2003). Additionally, the right panel of Figure 7 presents the shocked H2 emission in 
IRAS 19410+2336 (also taken from Beuther et al. 2003). We find many different CS emission 
features throughout the regions (CSl to CS9 in Fig. 7). However, the overall morphology 
differs significantly from the single-dish dust continuum data. The CS emission peaks fall 
into two classes, those in the near vicinity of the mm peaks (triangles in Figure 7; CSl, CS6 
and CS7), and the rest of the emission peaks which are all offset from the main mm cores. 
Comparing these offset CS peaks with the known outffows in the region, it appears that they 
are all associated with one or other outffow feature: CS2 is at the wake of outffow (B), CSS 
and CS4 appear to be associated with the outffows (G) and (H), CSS is spatially associated 
with outffow (J) and also with a strong H2 feature. Finally, CSS and CS9 both seem to be 
associated with outffow (F). While CSS is again at the wake of the outffow, CS9 is within 
the outffow lobe where there is only very weak CO outffow emission (Beuther et al. 2003). 
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Zooming into the main peak regions in the south (Fig. 8) and in the north (Fig. 9, left) 
with the higher resolution of the CS(5-4) data, we find that, although there is emission in the 
near vicinity of the mm peaks, the CS emission is offset from the main peaks. As pointed out 
in §3.2, the positional accuracy between the various maps is extremely reliable because of the 
simultaneous interferometric observations of the CS line and the mm continuum emission. 
While we find emission at all velocities in the north, the southern CS(5-4) emission is mainly 
found at blue velocities and the velocity of rest {Av [14,22] km/s). There is barely any red 
CS(5-4) emission in the southern core at all (Fig. 6). Comparing it again with the outflow 
emission (Beuther et al. 2003), we see that the blue outflow wing (C) follows the same 
direction as the blue CS emission. The CS spectra observed toward CSl (Fig. 6) somewhat 
resemble the infall spectra observed in CS toward the massive star-forming region W51 (Fig. 
8 in Zhang et al. 1998). However, as our data suffer strongly from missing flux problems and 
the negative dip at 23 km/s is observed toward positions offset from the peak CSl as well, 
we attribute this feature to observational artifacts due to missing short spacings and refrain 
from a physical interpretation. 

The image is less clear in the north, the CS peaks are offset from the mm peak emission 
(Fig. 9), and three outflows emanate from that region. Nevertheless we cannot unambigu- 
ously attribute the CS emission to any of these three outflows. So, it is very possible that 
other reasons contribute to the CS-mm offsets in the north as well, e.g., optical depth and/or 
freeze out of CS onto grains in cold core centers (see, e.g., Bergin & Langer 1997; Tafalla 
et al. 2004). In order to better understand the differences between the CS(5-4) and (2-1) 
transitions, we confined both datasets to the same uv-range and re-imaged these data. The 
missing flux problems are again severe (see §3.2) but we could produce integrated CS maps 
of both transitions (Fig. 9, right). Both transitions peak offset from the continuum, but it 
appears that the (2-1) transition is even further offset from the mm continuum peak than 
the (5-4) transitions. Different explanations for these observations arc possible: it could be 
an opacity effect with CS(5-4) turning optically thick closer to the center than the CS(2-1) 
observations. However, physical reasons like higher densities and/or temperatures toward 
the core-center can produce this observational effect as well. Unfortunately, based on this 
data we cannot differentiate between the different scenarios - independent observations of 
temperature and density tracers as well as of optically thin isotopomers are necessary to get 
a hand on the different processes taking place in this region. 

Similarly as in §3.2, we refrain from temperature and column density estimates based 
on the CS data because the missing flux deteriorates every such an estimate too strongly. 
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4. Discussion 

4.1. Line emission 

The search for and investigation of disks in massive star formation still remains one of the 
missing links for a more detailed understanding of high- mass star formation. Theories predict 
disks in massive star formation (e.g., Yorke & Sonnhalter 2002), and other observations claim 
indirectly that disks need to exist because more and more collimatcd massive outflows have 
been found (e.g., Bcuthcr ct al. 2002a). However, direct disk evidence in high-mass star 
formation is very scarce (Ccsaroni et al. 1997; Zhang et al. 2002) and yet not completely 
unambiguous. Therefore, it is necessary to establish a procedure to study the dense central 
cores and potential disks in massive star formation. We believe that thermal molecular line 
emission from dense gas tracers which are not prominent in the ambient medium is the most 
promising approach. 

One of the results of this paper is that CS and CN are not well suited for disk and dense 
gas studies in these sources. In spite of being a rather negative result, it needs to be stressed 
and we recommend the use of different molecules for future massive disk studies. 

High optical depth probably plays a role for CS avoiding the dense centers of the regions. 
It remains to be tested whether rarer isotopomers like C^'^S or ^^CS might trace the dense core 
centers and potential disks in a better way. In addition, the observations demonstrate that 
the CS emission is dominated to a large degree by the various outflows in IRAS 19410+2336. 
This CS emission is still associated with the large scale dust emission observed with single- 
dish instruments which likely reflects the similar range of gas densities needed to produce 
dust continuum and CS line emission. 

The situation in slightly different regarding the CN emission in IRAS 20393+3952: it 
also avoids the core centers but optical depth effects appear to be negligible. UV radiation is 
necessary to produce significant CN emission (e.g., Simon 1997), and CN is detected toward 
low-mass disks due to photo-dissociation by the star on the surface of the disk (Dutrey et al. 
1997; Aikawa & Herbst 2001). Two main reasons are likely to cause the lack of CN emission 
toward IRAS 20293+3952: (a) the source is too deeply embedded to excite the surface layers 
of the potential disk, or (b) the source is in such an early evolutionary stage that it does 
not produce enough photons. With the data available so far, we cannot properly distinguish 
between the two possibilities. 

These observations outline the difficulties of studying potential disks in massive star 
formation. We have to look for other molecular transitions to properly investigate massive 
disks. While NH3 sometimes does a good job (Zhang et al. 2002) it is often confused by 
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the ambient gas. The sulphur-bearing species SO was observed toward the central driving 
source of HH 80-81 (Gomez et al. 2003) but it is also found toward the outflows in massive 
star-forming regions (e.g., IRAS 18089-1732, Beuther et al. 2004c) implying that SO can be 
ambiguous as well. The most promising approach might be to focus on real hot-core-tracing 
molecules like CH3CN in IRAS 20126+4104 (Cesaroni et al. 1997, 1999), or HCOOCH3 
observed toward IRAS 18089—1732 with the Submillimeter Array (SMA) (Beuther et al. 
2004a). Prom an observational point of view, the SMA with a broad bandwidth of 2 GHz 
in two sidebands (4 GHz of spectral data) offers a unique opportunity to observe many 
interesting lines simultaneously and thus compare different molecular transitions within the 
same observation (Ho et al. 2004). This way, we should soon get a better census on the 
disk-tracing molecules in massive star formation. 



4.2. Continuum emission 

The dust opacity indices (5 we derive from the imaging as well as from uv-fits are 
significantly lower than the canonical value of 2 (Hildebrand 1983). These results resemble 
similar recent findings of low (3 in other massive star-forming regions (Beuther et al. 2004a; 
Kumar et al. 2003). However, due to insufficient signal-to-noise ratio and uv-coverage, the 
previous study toward IRAS 18089—1732 was only capable of determining decreasing (3 with 
increasing baseline length (Beuther et al. 2004a). Thus, their data could not tell whether 
the decreasing {3 is due to unresolved clumpy structure throughout the core or an actual 
decrease in (3 toward the core center. The new data toward IRAS 20293-1-3952 now resolve 
different spectral indices for the outer regions as well as the core center. It is not yet clear 
whether we simply have two different components - a dust core with constant (3 plus either 
an optically thick disk or an unresolved hypercompact Hll region - or a gradient of varying 
13 from the outside to the core center. Maybe both scenarios contribute to the observed 
spectral index variation. A gradient in (3 would be similar to results that have been reported 
in low-mass systems by Hogerheijde & Sandell (2000). Such a decrease of (3 could either be 
due to increasing opacity toward the core center or grain growth within circum-protostellar 
disks (see also Beckwith et al. 2000). Future work in that field would be required to cover an 
even broader range of spatial scales, and it will be interesting to determine at what distance 
from the core center (3 reaches the canonical value of 2, if it reaches it at all. 

The finding that a is lower toward the outflow than toward the core center is different 
from the recent results from Gueth et al. (2003) in the low-mass outflow LI 157, who reported 
an average spectral index of a ~ 5 toward the outflow. Their analysis is based on bolometer 
data over very broad band-passes (~ 80 GHz) and thus line contamination can introduce 
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artificial flux increases in one or both bands. But even taking such artifacts into account, 
Gueth et al. (2003) conclude that in L1157 a is higher toward the outflow than toward 
the core. Their observations point toward smaller dust grains and/or a different chemistry, 
which is consistent with models predicting the evaporation of icy dust grains during shock 
interactions of the outflow with the surrounding dust. 

These models obviously fail toward mmla. However, thermal ionized wind models 
predict spectral indices between 2 and —0.1 (Reynolds 1986). A typical spherical wind 
results in a ~ 0.6 whereas highly collimated jets have even lower values around 0.25, and in 
the optically thin case —0.1 (Reynolds 1986). Therefore, the derived «(mmla) ~ 0.8 within 
the outflow region of IRAS 20393+3952 is consistent with standard models of collimated 
ionized winds. 

We very much appreciate the comments from an anonymous referee which helped to 
improve the manuscript. We would like to thank the staff of IRAM in Grenoble for help 
during the observations and the data reduction process. H.B. acknowledges financial sup- 
port by the Emmy-Noether-Program of the Deutsche Forschungsgemeinschaft (DFG, grant 
BE2578/1). 
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Table 1. Observation parameters 





IRAS 20293+3952 


IRAS 19410+2336 


Date 


winter 2002/2003 


winter 2001/2002 & 2003/2004 


R.A.[J2000] 


20:31:12.9 


19:43:11.4 


Dec.[J2000] 


+40:03:22.9 


+23:44:06 


^'Isr 


6.3km/s 


22.4 km/s 


Distance" [kpc] 


2.0 


2.0 


Config. 


AB 


BCD 


Lines 


CN N=l-0 & N=2-l 


CS J=5-4 & J=2-l 


Freq.[GHz] 


113.385 & 226.656 


97.981 & 244.936 


Freq. res. 


0.312 MHz 


0.039 MHz 


Velo. res. 


0.83 & 0.41 km/s 


0.12 & 0.05 km/s 


^(1 mm) 


0.76" X 0.61"(RA.37°) 


1.52" X 1.03"(P.A.63°) 


^(3 mm) 


1.46" X 1.21"(P.A.32°) 


5.31" X 3.42"(P.A.77°) 


Mosaic offsets ["] 




-2/-2 -12/-8 9/2 -1.5/9.5 






-9.5/55 -19.5/49 -9.5/43 
-19.5/37-19.5/25 -10/17 


Phase calib. 


2013+370, 2023+336 


2023+336, 1923+210 


Amp. calib. 


MWC349 


MWC349 


3(Trms(1.3mm) 


3.6 mJy 


9mJy'' 


3crrms(2.6mm) 


1.2 mJy 


0.75 mJy'' 



^The kinematic distance ambiguities are solved by associating the regions via the near- 
and mid-infrared surveys 2MASS and MSX on larger scales with sources of known distance 
(Bontemps, priv. comm.). 

^Data presented in Beuther & Schilke (2004) 
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Table 2. Continuum data in IRAS 20293+3952 



Source 


Offset 

n 




Cint 

^2.6mm 

[mJy] 


Qpcak 
2.6mm 

[mJy/beam] 


cint 
1.3mm 

[mJy] 


ripeak 
1.3mm 

[mJy/beam] 


a 


-^^3mm 

[Mo] 


^Imm 

[Mo] 


mml 




0/0 


28.0 


16.1 


118.0 


42.3 


< 3.4'' 


2.8 


1.6 


mmla 


-2.7/ - 


- 1.3 


2.9" 


2.9 


4.8 


4.8 


~ 0.8 






mm2 


-10.2/ - 


10.7 


8.6 


5.6 








0.9 




mm3 


-5.6/ - 


12.4 


3.6" 


3.6 








0.4 





Source is unresolved. 
'See the main text. 
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Fig. 1. — IRAS 20293+3952: the grey-scale shows the 1.2 mm large-scale emis- 
sion obtained with the MAMBO-array (Beuther et al. 2002b), the grey-scale levels 
are 18 to 336mJy/(ll"beam) in 36mJy/(ll"beam) steps (18(36)336 mJy/(ll"beam)). 
The contours present the PdBI observations: white contours the 2.6 mm continuum 
(1.2(1.2)15.6 mJy/beam), and black contours the integrated CN(l-O) emission [Av 
[4,12] km/s, 29(42)407 mJy/beam). The arrows sketch the four outflows previously observed 
by Beuther et al. (2004c), and the beam is shown at the bottom left. The small box outlines 
the region presented in Figure 4. 

Fig. 2. — Millimeter Continuum emission in IRAS 20293-1-3952. The left and middle panels 
show the 2.6 mm and 1.3 mm emission, respectively. The contouring is done in 3a steps 
(1.2 mJy and 3.6 mJy, respectively). The right panel presents in grey-scale the spectral 
index map between the two bands and in contours the 1.3 mm continuum emission. For the 
right panel, the data of both bands were selected to cover the same uv-range (in kA). 

Fig. 3. — The 1.3 and 2.6 mm continuum data in the uv-plane and the corresponding two- 
component flts. The top data-points (triangles) are the 1.3 mm data and the bottom data- 
points are the measurements at 3mm (squares). The dotted lines show the fitted point 
source and the solid lines the combined fits consisting of a Gaussian plus a point source. 
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Fig. 4. — IRAS 20293+3952: the grey-scale with dotted contours shows the 1.3 mm con- 
tinuum emission (3.6(3.6)39.6 mJy/beam), and the hned contours present the integrated 
CN(2-1) data {Av [4,9]km/s, 111(74)703 mJy/beam). The beam is shown at the bottom 
right. 

Fig. 5.— CN(l-O) spectrum in IRAS 20293+3952 integrated over a region of AR.A.[-5, 1]" 
and ADec.[— 5, 1]", the two difTerent transitions are labeled at the top. 

Fig. 6.— IRAS 19410+2336: CS(2-1) and CS(5-4) spectra in full and dashed lines, respec- 
tively. The spectra are taken toward the positions CSl (within 6 arc-sec^) and CS6 (within 
9 arc-sec^) (Fig. 7). The CS(5-4) spectrum toward CSl is scaled up by a factor 2. 

Fig. 7. — IRAS 19410+2336: Left panel: the grey-scale shows the 1.2 mm large-scale 
emission obtained with the MAMBO-array (Beuther et al. 2002b), the grey-scale levels are 
85(85)765 mJy/ (ll"beam). The contours present the integrated CS(2-1) PdBI emission {Av 
[17,27] km/s, 0.2(0.2)1.8 Jy/beam). The arrows, eUipses, triangles and labels mark the out- 
flows and mm sources as outlined in Beuther et al. (2003). The 3 mm beam is shown at the 
bottom right. Right panel: the grey-scale now presents the infrared shocked H2 emission 
taken from Beuther et al. (2003). The contours again show the integrated CS(2-1) emission, 
and the numbers label the different CS emission features discussed in the main body. 

Fig. 8. — IRAS 19410+2336: the grey-scale shows the 1.3 mm continuum data already pre- 
sented in (Beuther et al. 2004c), the grey-scale levels are 7(7)63 mJy/beam. The contours 
outline the blue CS(5-4) emission in that region {Av [14,19] km/s, 0.5(0.5)4.5 Jy/bcam). 
The arrows sketch outflow C from Beuther et al. (2003), the synthesized beam is shown at 
the bottom right. 

Fig. 9. — IRAS 19410+2336: (left): the grey-scale shows the 1.3mm continuum data 
already presented in Beuther et al. (2004c), the grey-scale levels are 9(3)50 mJy/beam. 
The contours outline the integrated CS(5-4) emission in that region {Av [17,24] km/s, 
0.7(0.7)6.3 Jy/beam). In the right image, we present the CS and 1.3mm continuum data 
imaged all with the same uv-range for better comparisons. The grey-scale shows the con- 
tinuum, full contours CS(5-4) and dashed contours (CS(2-1), the contouring for all three is 
from 5 to 95% (step 10%) of their peak emissions, respectively. The synthesized beams are 
drawn at the bottom left of each panel. 
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